Lattice versus 2PI: 2d O(N) model at nonzero T* 
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The thermodynamics of the O(N) model in 1+1 dimensions is studied 
applying the CJT formalism and the auxiliary field method as well as fully 
nonperturbative finite temperature lattice simulations. The numerical re- 
sults for the renormalized mass of the scalar particles, the pressure and the 
trace anomaly are presented and compared with the results from lattice 
simulation of the model. We find that when going to the two loop order we 
observe a good correspondence between the CJT formalism and the lattice 
study. 
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1. Introduction 

The two dimensional O(N) nonlinear sigma model has many interest- 
ing features in common with four-dimensional non-Abelian gauge theories, 
which makes it useful to study as a toy model for QCD [IJ [2j O HI [5]. 
For instance, the coupling constant is dimensionless, therefore the theory is 
renormalizable. Besides, this model is asymptotically free and has a dynam- 
ically generated mass gap. Another interesting property is the conformal 
invariance: In two dimensions the nonlinear sigma model is classically scale 
invariant. However, on the quantum level a scale is introduced due to renor- 
malization of the quantum corrections. Furthermore, for N = 3 the model 
exhibits instanton solutions. We start with the usual expression for the 
generating functional at nonzero temperature 
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where g is the coupling constant and jCq is a free Lagrangian 
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The fields are restrained by the condition <I> 2 = iV/g 2 which is incorporated 
by the delta function. The nonlinear constraint enforces the thermodynam- 
ics of the model on an N — 1 dimensional hypersphere and induces the 
interactions between the fields. Using the mathematically well-defined (i.e., 
convergent) form of the usual representation of the functional ^-function 
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the generating functional and the corresponding Lagrangian of the O(N) 
nonlinear model can be rewritten as follows [6j 
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where a is an auxiliary field serving as a Lagrange multiplier. 



2. Analytic calculations 

In order to study the thermodynamic behavior we apply the CJT for- 
malism. Within this formalism we obtain the following expressions for the 
renormalized effective potential and for the renormalized gap equation to 
one-loop order 
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where fi is the renormalization parameter, m is the vacuum mass, M is the 
temperature dependent mass and g 2 en is the renormalized coupling constant. 
The details of the computation can be found in [7]. Next, we can calculate 
the thermodynamic pressure which is, up to a minus sign, identical to the 
minimum of the effective potential 
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Fig. 1. The analytic result for the mass of the scalar field as a function of T/m. 



As was shown in [5], to two-loop order one must apply numerical methods 
in order to regularize the effective potential at finite temperature. The final 
result is rather lengthy and is given in [5] as well as in [7]. 



3. Lattice simulation 

In this section we summarize the thermodynamic approach applied by 
finite temperature lattice simulations. The Euclidean, discretised action 
takes the form of a Heisenberg model, 

s = pY,0--%-%) > (3) 

(iJ) 

where the sum runs over all bonds of a 2-dimensional lattice, s*j are 3- 
dimensional unit vectors in internal space and /3 = N/g 2 . The corresponding 
partition function reads 




The system at finite temperature T is realized by making the time-like 
extent of the lattice finite and consisting of Nt sites, with periodic boundary 
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conditions in that direction; denoting with a the lattice spacing, we have 
aN t = l/T. 

In order to evaluate the pressure we use the "integral method" [8] and 
additive renormalization 



where £ e = s$ • Sj+g. This requires the knowledge of the "beta function" 
d/3/d In T = —ad/3 /da, which is extracted nonperturbatively, measuring the 
running of the coupling with the scale and using a suitable parametrization 
of the data. Note that the system at T = is numerically approximated 
with a large enough square system, that is N x = Nt- 

The running of the scale a(/3) is needed to convert the temperature 
in physical units. This is realized by computing the spin-spin correla- 
tion function in a zero-temperature system over a wide range of couplings 
C{an) = {si ■ Si+ n e), from which we obtain the zero-temperature mass m. 
Distances on the lattice are in units of a, therefore the masses we measure 
are the adimensional quantities am; once we fix the corresponding physical 
value we have the function a{j3) in units of 1/m. 



In this subsection we compare the analytic results with lattice simula- 
tions. We choose N = 3, corresponding to a system of three scalar fields. 
The vacuum mass of the scalar field m is the only dimensionful parameter 
of the model. Therefore we plot all thermodynamic observables in units of 
m. 

Figure [T] shows the analytic result for the temperature dependent renor- 
malized mass of the scalar fields. An important observation is that it has 
a similar behavior as the gluon mass in the deconfined phase, e.g. ref. [U] 
and refs. therein. The model exhibits dimensional transmutation just like 
QCD, meaning that at zero temperature there is a nonvanishing mass gap, 
which is generated due to renormalization of quantum corrections. Besides, 
at high T the temperature dependence of the mass can be approximately 
parametrized by T/logT. This is exactly the same behavior that one ob- 
serves for the gluon mass in the deconfined phase. 

Figure [2] shows the thermodynamic pressure. The uncertainties are very 
small and completely overcome by systematics. The high temperature limit 
is better described in the two-loop calculations, whereas at low temperatures 
there is a better correspondence between the lattice simulations and the 
one-loop calculations. We find that at very high temperatures one degree 
of freedom becomes effectively removed, and the pressure approaches the 





4. Results and Discussion 
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Fig. 2. The pressure as a function of T/m in one-loop (blue) and in two-loop (red) 
approximation compared to lattice simulations for different values of (3 = N / ' g 2 . 



limit of a non-interacting gas of iV — 1 bosons, P/N = (N — l)irT 2 /6N, 
which is indicated by the straight dashed line in the upper right corner. 
This is an immediate consequence of the asymptotic freedom and of the 
nonlinear constraint. One can understand this behavior by remembering 
that we start with a free Lagrangian for N scalar fields. Introducing the 
nonlinear condition, <I> 2 = N/g 2 , the thermodynamics is constraint on a 
two dimensional sphere which effectively removes one degree of freedom. 
Analytically this removal of one degree of freedom can only be described in 
the 2-loop approximation. 
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